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ABSTRACT 

We suggest that planets, brown dwarfs, and even low mass stars can be formed by fragmentation 
of protoplanetary disks around very massive stars (M > lOOM©). We discuss how fragmentation 
conditions make the formation of very massive planetary systems around very massive stars favorable. 
Such planetary systems are likely to be composed of brown dwarfs and low mass stars of ~ 0.1 — 
0.3M Q , at orbital separations of ~ few x 100 — 10 4 AU. In particular, scaling from solar-like stars 
suggests that hundreds of Mercury-like planets might orbit very massive stars at ~ 10 3 AU, where 
conditions might favor liquid water. Such fragmentation objects can be excellent targets for the 
James Webb Space Telescope and other large telescopes working in the IR bands. We predict that 
deep observations of very massive stars would reveal these fragmentation objects, orbiting in the 
same orbital plane in cases where there are more than one object. 

Subject headings: protoplanetary disks — planets and satellites: formation — stars: formation — stars: 
massive — stars: low-mass — (stars:) brown dwarfs 



1. INTRODUCTION 

Planets were found around objects which had been 
considered not to be able to support planets, such as 
around the pulsars PSR182910 (Bailes et al. 1991) and 
PSR1257+12 (Wolszczan & Frail 1992), where they are 
believed to be formed after the supernova explosion. In 
other cases planets can survive the post main sequence 
evolution of their parent stars. Examples include a plan- 
etary mass companion in orbit around V391 Pegasi (Sil- 
votti et al. 2007), and the close planet around HD 149382 
(Geier et al. 2009), both are extreme horizontal branch 
stars. These planets survived the red giant branch phase 
of their parent star. There was also a tentative detec- 
tion of a planet around a white dwarf (e.g., Mullally et 
al. 2008). Theoretical models even predict formation of 
second generation planets in disks around MS stars ac- 
creting mass from their AGB companions (Perets 2010). 
In this paper we are studying the possibility of finding 
planets, brown dwarfs (BDs) and low mass stars orbit- 
ing in the same plane around very massive stars (VMS, 
M» > 100M Q ). 

For planets to exist around VMSs their orbital dis- 
tance must be very large. A relevant and famous extra- 
solar planet was observed around the ~ 2M© star Foma- 
lhaut (Kalas et al. 2008). The orbital distance of the 
observed planet Fomalhaut b is ~ 115 AU, and its mass 
was constrained to be < 3Mj (Chiang et al. 2009). Fo- 
malhaut b teaches us that planets can be formed at large 
distances from their parent star. 

According to common view, gas giant planets can be 
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formed through two main channels (e.g., Boley 2009; 
Dodson- Robinson et al. 2009): accretion of planetesimals 
into a core, followed by gas accretion (e.g., Pollack et al. 
1996; Kenyon & Bromley 2009; Brauer et al. 2008), and 
instabilities in the protoplanetary disk (PPD) that causes 
fragmentation (e.g., Cameron 1978; Mayer et al. 2002; 
Boss 1997, 2006; Brauer et al. 2008). Dodson-Robinson 
et al. (2009) showed that giant planets in wide orbits of 
r > 35 AU form by instabilities in the PPD. Stamatel- 
los & Whitworth (2008) found that fragmentation does 
not occur below r < 40 AU, a result in agreement with 
Dodson-Robinson et al. (2009), considering the some- 
what different parameters used. 

Rice et al. (2003) found using hydrodynamic simula- 
tions and N-body orbit integrations that when a O.IMq 
PPD around a 1M star fragments, the most massive ob- 
ject to be formed is a 0.0075Mq BD. As the simulation 
of Rice et al. (2003) is scale-free, it immediately implies 
that a PPD around a IOOMq star is expected to fragment 
into a low mass star of up to 0.75M©. A more detailed 
estimate carried in section[2]gives a somewhat lower mass. 

In this paper we examine the possibility that plan- 
ets, BDs, and in particular low mass stars can be formed 
around a VMS. We start with two basic assumptions: (1) 
Very massive single and binary stars posses circumstellar 
and circumbinary disks similar to those of low mass stars. 
(2) We can scale the constraints and conditions from low 
mass stars to VMSs. This assumption is based in part 
on the results of Kratter & Matzner (2006) who studied 
massive PPDs. Although they did not study explicitly 
the regime explored here, some of their parameters over- 
lap with those studied by us. Following these basic as- 
sumptions, in section [2] we derive scaling relations, and 
show that fragmentation around VMSs can occur, and 



produce fragmentation objects as massive as low mass 
stars. In section [3] we discuss the observational signature 
of the proposed fragmentation objects. We summarize 
our results in section [5] 

2. DISK FRAGMENTATION AROUND A 
VERY MASSIVE STAR 

We wish to determine the properties of a fragmented 
protoplanetary disk (PPD) around a VMS; we scale by 
M„ = 100M Q . The PPD is assumed to be vertically 
optically thick (i.e., perpendicular to the disk plane) at 
the distance where it fragments (Kratter et al. 2010). 
The surface density profile of the PPD is assumed to have 
the same form as in low mass stars 



E(r) = Eo (liu) P gcm " 



where Eo and p are constants (Weidenschilling 1977; 
Hayashi 1981; Nero & Bjorkman 2009). The power p 
varies from one PPD to the other in the range p — 0.5 — 3 
(e.g., Nero & Bjorkman 2009), but more commonly in the 
range p = 1 — 3/2. For early stages of PPD evolution, 
where steady accretion from the PPD to the star takes 
place, an analytical solution to the evolutionary equation 
of the PPD (e.g., Frank, King & Raine 2002) gives p = 1. 
At later stages the accretion to the star decreases and 
the value of p increases, until accretion stops and then 
an analytical solution to the evolutionary equation of the 
PPD gives p — 3/2. There are hints that practically, for 
more massive stars a smaller value of p should be used. 
For the solar protoplanetary nebula a value of p = 3/2 
is usually used (Weidenschilling 1977; Hayashi 1981), but 
other values, such as p = 2.168 (Desch 2007), were also 
suggested. The models of Vaidya et al. (2009) for PPDs 
around 10 — 37Mq stars were best fitted with p — 1.1. 
The stellar disk structure analysis of Andrews & Williams 
(2007) based on high resolution submillimeter continuum 
survey of circumstellar disks in the Taurus- Auriga and 
Ophiuchus-Scorpius star formation regions, also supports 
a value of p ~ 1. For the solar PPD, Eo is taken to be in 
the range of 1700 g cm" 2 (Hayashi 1981) to 4200 g cm" 2 
(Weidenschilling 1977). For the disk around Fomalhaut 
the model of Nero & Bjorkman (2009) implies a value of 
E ~ 10 3 g cm -2 . 

Our goal is to determine the range of the orbital dis- 
tances r v where fragmentation around VMSs occurs, and 
the outer radius of the PPD, Rd. For that we start with 
results from studies of PPDs around low mass stars, and 
our assumptions of similarity to VMS. We impose four 
conditions for fragmentation objects to be formed around 
VMSs: 

(1) Planets are formed by fragmentation in the PPD 
where the surface density is ~ 10 — 100 g cm -2 (e.g., 
Nero & Bjorkman 2009). We assume that fragmentation 
in PPDs around VMSs occurs at the same surface density. 

(2) The solar PPD mass was ~ 1 — 20 per cent of the 
mass of the Sun (Weidenschilling 1977; Hayashi 1981; 
Kuchner 2004; Desch 2007). Again, we take approxi- 
mately the same fraction to hold in our study, which 
translates to a total PPD mass of ~ 1 — 20Mq. This 



mass leads to a relation between the surface density and 
the outer radius of the PPD, Rd 



E(r)27rr dr, 



(2) 



where we assume that the inner radius of the PPD is 
small and can be neglected (for p < 2 in equation ifljl). 

(3) In order for the PPD to fragment, two conditions 
have to be satisfied (e.g., Kratter et al. 2010). (3.1) 
Toomre's parameter (Toomre 1964) should satisfy 
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where Q is the Keplerian angular velocity, c a — \JkT j fj,m p 
is the speed of sound, m p is the proton mass, \x is the 
molecular weight, and k is the Boltzmann constant. (3.2) 
The radiative cooling time (in the optically thick PPD) 
must be shorter than the orbital time (Gammie 2001), 



^cool ^ 



(4) 



so the fragment can cool in less than about one orbital 
period. According to the model of Kratter et al. (2010) 
there is a minimal orbital distance which satisfies both 
conditions for fragmentation , the fragmentation distance. 
For a 1.5M0 star, Kratter et al. (2010) found this radius 
to be ~ 70 AU. An approximatively similar fragmenta- 
tion distance was obtained by Clarke (2009). Based on 
the results of Kratter et al. (2010) and our assumptions, 
we scale with a VMS to derive the fragmentation distance 
appropriate for our studied parameter space 



r f ~ 280 



Al, 



100M 



AU. 



(5) 



The same fragmentation distance can also be determined 
from the condition that the vertical optical depth is r ~ 1 
(Kratter et al. 2010). As can be seen the fragmentation 
distance weakly depends on the mass of the central star. 
As for VMSs the mass-luminosity relation is A/* oc L», 
where /3 ~ 0.3, so the fragmentation distance very weakly 

depends on the luminosity, r/ oc L* . 

(4) Observations show that planets form quite close 
to the fragmentation distance. For example, Fomalhaut 
b, the outer planets in the triple-planetary system HR 
8799, and the potential protoplanet associated with HL 
Tau, were all formed outer to, but quite close to the frag- 
mentation radius, at distances of ~ (1 — 3)r/. For our 
studied cases we approximate this range as ~ (1 — 5)r/ 
which translates to ~ 300-2000 AU if we consider VMSs 
with masses up to ~ 3OOM0. We would therefore expect 
to find the inner fragmentation objects approximately 
within this range. Outer fragmentation objects can re- 
side at larger orbital distances, up to ~ 3000 AU. 

We build three calibrations which obey all the con- 
ditions listed above, with the main difference being the 
power p of the surface density profile given in equation 
0- First, we examine the implications of the demand 
that for M, — 100M© the fragmentation object is formed 
at r v = 1000 AU and that the surface density there is the 



minimal value allowed by condition (1), E p = 10 g cm 
We therefore obtain our three density profiles 
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As can be seen from equation (|6]) and the four condi- 
tions, there is some freedom in the values of the different 
numerical factors, and the set of numbers is not unique. 
However, the possible range of numerical factors do not 
change much our conclusions. The ratio between the or- 
bital distances r p where fragmentation around VMSs oc- 
curs and the fragmentation distance 77 is taken to be 
constant \ = r p /rf. This constant ratio is set such that 
for M» = 1OOM we shall get r p = 1000 AU. As can be 
seen from equation (JS| this ratio is x — 3.5. We therefore 
get, for all three calibrations 
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(7) 



Note that according to our conditions (3) and (4) the 
fragmentation objects can be formed at larger distanced 
than r p = 1000 AU, up to 



1500 



AU. 



(8) 



V 100M Q 

The PPD is assumed to posses a total mass of md = 
0.2M, , the maximum allowed by condition (2) , and there- 
fore extends up to 



Rd 




■p=l 



;p=S 
iP= I- 



(9) 



The most simple approximation for the mass of the 
fragmentation object, using the models of Goodman & 
Tan (2004) and Kratter & Matzner (2006) gives 
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(10) 

where E p = E(r p ) was calibrated to a range compatible 
with planet formation as discussed above, fl p = fl(r p ), 
Q p = Q(r p ) and c SjP = c s (r p ). In the second equality of 
equation (|10p we used equations ©, (J7J), and the expres- 
sion for the speed of sound. 

We will also check an approximation from another 
model. As the PPD is assumed to be optically thick at the 
distance where it fragments, the mass of the fragmenta- 
tion object can be approximated as in Nero & Bjorkman 
(2009) 
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Note that in equations (|10p and (|1 If) E p and r p are 
related by equation © and cannot be individually cal- 
ibrated. Substituting equations ([6| and (0 in equation 
(fTTj) we get 
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namely the dependence of m / in the mass of the VMS is 
the same for our three calibrations. Note that the typical 
mass of the fragmentation object is on the boundary of 
BDs and low mass stars. 

Let us estimate the temperature at the radius r p where 
fragmentation takes place. The simplest estimate is to use 
Toomre's parameter (equation ©), and to substitute the 
explicit expression for the sound speed. This gives for the 
temperature 



T(r p ) = 
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(13) 



K. 



where in the second equality we used equation ([7]) to elim- 
inate the dependance on the stellar mass (which comes 
from Q p ). 

We will make another estimate of the temperature ar 
r p . Using the Q-disk model for PPDs as in Dullemond 
et al. (2007), we can estimate the intrinsic temperature 
at the center of the PPD which depends on the radius as 
T c (r) oc r~5. The PPD is assumed to be geometrically 
thin, and the temperature is assumed to result from the 
internal viscosity of the PPD, and not from the radiation 
of the parent VMS. Namely, the PPD is horizontally op- 
tically thick. The model of Dullemond et al. (2007) uses, 
as in our second calibration, a surface density profile of 
E oc r , explicitly written as 



E(r) 
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(14) 



where M is the mass accretion rate of the PPD and a is 
the a-disk parameter (Shakura & Sunyaev 1973). Vaidya 
et al. (2009) found that a ~ 0.1 is an optimal value in 
order for a PPD to be on one hand stable against com- 
plete fragmentation due to thermal effects caused by high 
viscosity (occurs for a > 1), and on the other hand not 
completely subjected to gravitational instability (occurs 
for a < 0.01). Namely, setting a ~ 0.1 well describes 
the situation observed in PPDs, where fragmentation oc- 
curs beyond some fragmentation distance, such as given 
in equation ([SJ. For detailed analytical calculations of a 
and the PPD's viscosity profile see Isella et al. (2009). 

Equating equation (|14[1 with the first calibration of 
equation ([S]), we find that the PPD's central (midplane) 
temperaturetemperature in the radius r p where fragmen- 



tat ion takes place is 
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(15) 

In the second equality we used equations (0 and the first 
calibration in equation © to eliminate the dependance 
on the E p , r p , and M„. 

We consider also the case were the central star radi- 
ation is not blocked but rather heat the disk at r v . An 
estimate from another model, by Dodson- Robinson et al. 
(2009), based on an older model of Chiang & Goldreich 
(1997) which takes the existence of dust into account, 
gives (when translated from stellar radius and effective 
temperature to luminosity) 

T(rg^ 200 (^ i (_^)-^_^_^ K , 

(16) 

where L* is luminosity of the VMS, and 9 is the flaring 
angle of the PPD at r — r p . 

The temperature estimates in equations (|13p and (|13p 
are compatible with estimates of temperatures where 
fragmentation usually takes place, and planets are formed 
< 100 K (e.g., Boss 1998, 2009; Dodson-Robinson et al. 
2009). The high temperature derived in equation (|16[) 
shows that if the radiation from the VMS is not blocked, 
then Q > 1, and there will be no fragmentation. 

Our most significant conclusion of this section is that 
in PPDs around VMS, low mass stars and BDs can be 
formed in the same way as planets are formed around 
solar-like stars. For example, we might find 2 or more 
low mass stars and BDs orbiting the parent star in the 
same plane as BDs and planets around low mass stars. 

Though in some cases our calculation can be refereed 
to as guidelines, and the exact numbers might be some- 
what different, we showed that using different approaches 
to calculate the physical properties of the PPD, we find 
that it is likely to fragment and form an object in the 
PPD. 

3. OBSERVATIONAL PREDICTIONS 

The ability to observe fragmentation objects depends 
strongly on their masses, and their orbital distance from 
the parent VMS. The fragmentation objects to parent 
star typical mass ratio is very small, ~ 10 -3 , and their 
typical orbital period is extremely long, > few x 1000 yrs. 
It is impossible to detect such objects by the Doppler shift 
method. If the fragmentation object is a planet it would 
only reflect the light of its parent star, while if it is a BD 
or a star it would produce its own luminosity, peaked in 
the infrared. 

In a relatively short time a VMS enters the Lumi- 
nous Blue Variable and Wolf-Rayet stages of evolution, 
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Fig. 1. — Dashed thin black line: theoretical blackbody flux 
(in units of 10 26 erg s — *A ) from a low mass star or brown 
dwarf with an effective temperature of 3000 K and a luminosity 
of 10 — 3 Lq in an early evolution stage. Solid thin blue line: 
theoretical model (Rice et al. 2010) for emission of a 3000 K 
low mass star (in units of 10 26 erg s -1 A 1 ). Solid thick red 
line: The flux ratio (xlO -7 ) between a (3000 K; 10 _3 L Q ) low 
mass star or brown dwarf and a (50000 K; 2 X 10 6 Lq) very 
massive star, both assumed to emit blackbody radiation. The 
James Webb Space Telescope (JWST) will observe in the band 
X = 0.6 — 26 (im, and can detect the fragmentation objects 
studied here. 

in which it expels a considerable fraction of its mass (e.g., 
Smith & Owocki 2006). Suppose it expels a fraction 
AM/A1, of its mass over a time scale longer than the 
orbital period. The orbital radius of the fragmentation 
object would move to a larger distance of 

— ™(i^)(^)(^r-< 

(17) 

where r p is the orbital distance where the fragmenta- 
tion object is formed (Eggleton 2006). If there are short 
phases of large mass loss the orbit will become eccen- 
tric, with the distance to apastron, where the companion 
spends most of the time, larger even. Low mass stellar 
companions might be photometrically detected at large 
orbital separations of ~ 2000— 10 4 AU, as the separation 
at a distance of, say, 2 kpc be ~ 1 — 5". For exam- 
ple, it might be possible to detect a low mass star with 
the James Webb Space Telescope (JWST). We take a 
low mass star (or a BD) with an effective temperature 
of 3000 K and luminosity of the order of ~ 1O _3 L0, and 
a VMS with an effective temperature of 5 x 10 4 K and 
luminosity of ~ 2 x 10 6 Lq, and find that the ratio be- 
tween the fluxes of the two stars at A = 1 /im, close to 
the wavelength where the fragmentation object peaks, is 
~ 10~ 7 . As the separation is large, this ratio should be 
detected by the JWST. This is depicted in Fig. Q] 

In the scenario described in section [2] more than one 
fragmentation object may be formed. If a few planets, 
BDs or low mass stars are formed, they are expected to 
(1) be located in the same orbital plane, where the PPD 
used to lie, (2) have circumstellar distances in the range 
of a ~ few x 100 — 10 4 AU, and (3) have a total mass of 
a few per cents of the parent VMS. We predict that ob- 



servations of VMSs with the JWST and other telescopes 
of high capabilities in the IR bands can reveal systems 
where these conditions are fulfilled. 

We note the dramatic difference in the evolutionary 
timescale of the parent star and the fragmentation ob- 
jects. A VMS is expected to go through its entire evo- 
lutionary path in only a few x 10 6 yrs (e.g., Maeder & 
Meynet 1994). But a low mass star would not even reach 
its zero-age main sequence by then (e.g., D'Antona & 
Mazzitelli 1994), neither a BD would start burning Deu- 
terium (e.g., D'Antona & Mazzitelli 1985). At the age of 
few x 10 6 yrs the temperature of developing BDs may be 
~ 2000 - 5000 K (e.g., Chabrier et al. 2000), approxi- 
mately in the same range as for low mass stars. 

The estimations of Kratter & Matzner (2006) suggest 
that the disk-born stars (with masses of ~ 1M© for M* = 
lOOM©) are formed near or somewhat outside r p ~ 100 — 
200 AU and therefore there is a very low chance to detect 
them. According to our estimate it is possible that the 
fragmentation objects reach up to ~ 10 4 AU, and thus 
we predict that these low mass stars and BDs can be 
relatively easily observed. 

Recent observations of Smith et al. (2009) reveal an 
intriguing potential candidate for a fragmentation object 
resembling the one we suggest. The source G in W49A 
appears to be a 45M© protostar surrounded by a disk 
from which it was formed. The PPD is extending up to 
a distance of ~ 6 x 10 4 AU, and the inner disc has been 
probably cleared out to a radius of ~ 1000 — 2000 AU. 
Within the disk, in the same plane, there appears to be a 
second star located ~4x 10 4 AU from the central source. 
Though these distances are larger by about one order of 
magnitude than our preliminary calibration, the common 
orbital plane may be an indication that this star formed in 
the disk, perhaps in fragmentation process as we suggest. 

Vlemmings et al. (2010) have discovered a strong mag- 
netic field around the massive protostar Cepheus A HW2, 
which also has a disk with a radius of ~ 600 AU surround- 
ing it. This discovery suggests that strong magnetic fields 
threading the PPD and maintaining a high accretion rate 
to the star, which is essential to its formation. It is not yet 
clear to us if these magnetic fields are helpful for planet 
formation as though they assist the formation of the star, 
they may as well moderate instabilities in the disk and 
make fragmentation more difficult. 

4. PHOTOEVAPORATION CONSIDERA- 
TIONS 

As we suggest that even low mass stars can be formed 
around VMSs, a point worth discussing is whether the 
radiation of low mass stars prevent their creation. Our 
answer to this question is most likely no. The track on 
the HR diagram in which low mass stars settle in the 
main sequence, starts in luminosity much higher than the 
luminosity the low mass stars have when they reach the 
main sequence. Thus, as the PPD fragments the low 
mass star has already passed the critical point in which 
its luminosity had the chance to prevent its creation. We 
therefore do not expect the radiation of the fragmentation 



object be an obstacle for the low mass star formation. 

Krumholz et al. (2009) performed 3D hydrodynamical 
simulations of the formation of a massive star via accre- 
tion from a disk. It was clearly shown that during the for- 
mation of a massive star, while it is still accreting through 
a PPD, the radiation it emits does not prevent accretion, 
and neither causes the PPD's destruction. Most of the ra- 
diation of the massive star is instead tunneled into polar 
directions, where it creates two low density polar bub- 
bles. Material from the cloud continues to be accreted 
onto the star, and when it comes from polar directions it 
flows along the walls of the bubbles and accreted onto the 
star. The PPD remains for a few x 10 4 yr, long enough 
for fragmentation to occur, and form a companion star in 
the disk. Kuiper et al. (2010) have recently developed a 
promising frequency dependent radiation transport simu- 
lation that is likely to strengthen the results of Krumholz 
et al. (2009). 

Radiation from massive star photo-evaporates their 
PPD, and lead to the formation of a disk wind. We will 
examine the degree to which photoevaporation by VMS 
radiation disrupt the PPD and prevents fragmentation. 
We use the results from Hosokawa & Kazuyuki (2009) 
who modeled the evolution of a newly formed accreting 
star reaching a final mass of IOOMq during its formation 
(their model MD3x3, presented in their figure 18). In 
this model of Hosokawa & Kazuyuki (2009) the mass ac- 
cretion rate is fixed at M acc = 3 x lO _3 M0yr _1 , and the 
VMS reaches its final mass in ~ 3.3 x 10 4 yr. From the 
model of Hosokawa & Kazuyuki (2009) we calculate all 
the basic parameters of the star as a function of its evolu- 
tion time. The first panel of Figure [2] shows the luminos- 
ity, mass, effective temperature and radius of the star as 
a function of its evolution time. The final luminosity of 
the star is T e = 58, 2007^, and its effective temperature 
is L — 2.41 x 1O 4 L<0. Assuming a blackbody radiation, 
we calculate the ionizing photon rate for this model (the 
results of Schaerer & de Koter 1997 show that this is an 
adequate approximation). At t = 1.7x 10 4 yr the ionizing 
photons rate rapidly raises, to 47riVii « 10 50 s _1 and then 
disk photoevaporation should be considered. 

We follow the model of Hollenbach et al. (1994) which 
related the ionizing photons rate to the mass loss rate 
expected from the PPD 

More detailed models (e.g. Font et al. 2004; Alexander et 
al. 2006) give similar values to those obtained by Hollen- 
bach et al. (1994). The second panel of Figure [2] shows 
the ionizing photons rate, the disk wind mass loss rate as 
well as some other quantities. 

To find the total mass lost from the PPD by pho- 
toevaporation we integrate the mass loss rate over the 
evolution time of ~ 3.3 x 10 4 yr and find 

= J M wind at ~ 2.5M Q . (19) 

This is much smaller than the mass of the PPD in our 
model, and can be regarded as a negligible loss of mass. 
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Fig. 2. — First panel: The luminosity, mass, effective tem- 
perature and radius of the star as a function of its evolution 
time, as inferred from the model of Hosokawa & Kazuyuki 
(2009). Second panel: Calculated quantities based on the 
model of Hosokawa & Kazuyuki (2009): Mass accretion rate, 
the Kelvin-Helmholtz time, accretion time from the disk, ion- 
izing photons rate, and mass loss rate by the disk wind (eq. 
I18jl . When integrating the mass loss rate by disk wind over 
the evolution time (equation H19H we find that only ~ 2.5Mq 
were lost from the PPD, and therefore it survived photoevap- 
oration. 



We therefore conclude that photoevaporation cannot de- 
stroy the disk during the relevant time for fragmentation. 



5. SUMMARY AND DISCUSSION 

Planets are found in a variety of orbital distances, 
around stars ranging in mass from low mass to massive, in 
both metal rich and metal poor environments, around sin- 
gle stars and binaries, and around stars in different stages 
of evolution (see the Extrasolar Planets Encyclopaedia, 
maintained by Schneider 2010, and references therein; for 
a recent statistical analysis see Mordasini et al. 2009). In 
the present paper we suggested that very massive stars 



(VMSs) of M, > 100M© are very likely to also harbor 
planetary systems. 

Under the assumption that fragmentation of proto- 
planetary disk (PPD) around solar like stars can be scaled 
to match the properties of PPD around VMS, we have 
found here that it is likely that in such PPDs, BDs and 
low mass stars of ~ 0.1 — O.3M0 are likely to be formed. 
We suggest that these fragmentation objects are orbit- 
ing the VMS in circumstellar distances in the range of 
~ few x 100 - 10 4 AU, on the same orbital plane. The 
low mass stars might be photometrically detected in the 
near infrared with modern ground telescopes using adap- 
tive optics in the IR, and modern space telescopes such 
as the JWST. 

As in low mass stars (e.g., Qian et al. 2010) it is 
possible that the massive fragmentation objects will be 
formed around binary systems of VMS in P-type orbits 
(outside the two components of the binary system). The 
fragmentation objects have very large orbital distances 
that ensure stability (e.g., Szebehely 1980; Dvorak et al. 
1989; Musielak et al. 2005). We note that simulations by 
Krumholz et al. (2009) have shown that fragmentation 
can be possible also if in the center of the PPD lie a mas- 
sive binary system. In the same simulation, the massive 
binary companion itself was formed from fragmentation 
objects in the PPD around the primary which merged 
into a massive stellar companion. 

Theoretical calculations based on the common stelar 
formation process limit the mass of the newly born star 
to M < 150M Q (e.g., Figer 2005). However, it seems 
that stars can reach higher masses, e.g., the Pistol star 
(Figer et al. 2004; Najarro 2005)) and 77 Car (Kashi & 
Soker 2010). These higher masses suggest that accretion 
from the PPD continues for a relatively longer time, with 
lots of mass processed through the PPD. Merger of a 
number of stars to form VMS also implies that there is 
lots of material in the surrounding of the newly formed 
star. The PPD itself is replenished by the surrounding 
gas cloud. Our study shows that such PPD are likely to 
leave fragmentation objects and debris (smaller objects 
and dust). 

Our calculations focused on large fragmentation ob- 
jects. Around solar-like stars these are planets, while 
around VMS these might be scaled to brown dwarfs and 
low mass stars. We can take another, somewhat specula- 
tive, step, and scale the Kuiper belt from the solar system 
to VMS. In VMSs, such objects might be ~ 100 times as 
massive as those in the solar system, and Kuiper objects 
are scaled to Mercury like planets. 

It is possible that such planets will sustain water in 
its liquid state. Crudely extrapolating the conditions for 
stable P-type orbits (e.g., Musielak et al. 2005) to a bi- 
nary mass ratio of M2/MI = 10 -3 implies that a planet 
can have a stable orbit at a > 1.6 times the binary sep- 
aration. Let us consider a VMS with a luminosity of 
L, ~ 2 - 5 x 1O 6 L0, for which a temperature similar to 
that on earth will be found at a circumstellar distance 
of a ~ 1500 — 2200 AU (depending on atmospheric con- 
ditions). Therefore, if the outer large fragment object 
(of mass ~ 0.1 - 0.2M Q ) is at a < 1500 AU, Mercury 



type planets might exist in an habitable zone, such a 
zone around VMS might contain hundreds of Mercuty 
like planets, and many more smaller objects. Collisions 
between these will form dust. Therefore, a zone of excess 
in IR emission might exist there. In a future paper we 
will examine ways to detect debris disks around VMSs. 

The PPDs around VMSs are expected to be huge and 
very massive relative to those around solar-like stars. 
These types of disks have favorable conditions for forma- 
tion circumstellar objects via fragmentation, rather than 
core accretion. Population III stars are thought to be of 
VMSs. If our results hold for very low metallicity popu- 
lations of VMSs, then it is not unreasonable to consider 
the formation of planets and brown dwarfs (and of course 
low mass stars) even around population III stars. Proba- 
bly not around the very first stars, but around those that 
contain a very low abundance of metals, i.e., late popula- 
tion III stars with metallicity of ~ 10 -6 times solar, and 
early population II stars. The VMSs explode eventually, 
and the fragmentation objects become unbound (if more 
than half the mass is lost in the explosion). We there- 
fore raise the possibility that the (rare) most metal-poor 
stars in the galaxy were formed in PPDs of the very first 
Population II VMSs. 

We thank Scott Kenyon, Kaitlin Kratter, Ken Rice, 
Nathan Smith, Ralph Neuhauser, Dimitris Stamatellos, 
and Bhargav Vaidya for helpful comments. This research 
was supported by the Asher Fund for Space Research at 
the Technion, and the Israel Science Foundation. 

REFERENCES 

Alexander R. D., Clarke C. J., Pringle J. E. 2006, MN- 
RAS, 369, 216 

Andrews, S. M. & Williams, J. P. 2007, ApJ, 659, 705 

Bailes, M., Lyne, A. G. & Shemar S. L. 1991 Nature 352, 
311 

Boley, A. C. 2009, ApJL, 695, 53 
Boss, A. P. 1997, Sci., 276, 1836 

Boss, A. P. 1998, Annu. Rev. Earth Planet. Sci., 26, 53 
Boss, A. P. 2006, ApJL, 637, 137 
Boss, A. P. 2009, ApJ, 694, 107 

Brauer, F., Dullemond, C. P., & Henning, T. 2008, A&A, 
480, 859 

Cameron, A. G. W. 1978, Moon Planets, 18, 5 

Chabrier, G., Baraffe, I., Allard, F., & Hauschildt, P. 
2000, ApJ, 542, 464 

Chiang, E., Kite, E., Kalas, P., Graham, J. R. & Clampin, 
M. 2009, ApJ, 693, 734 

Chiang, E., & Goldreich, P. 1997, ApJ, 490, 368 

Clarke, R. J. 2009, MNRAS, 396, 1066 

D'Antona, F. & Mazzitelli, I. 1985, ApJ, 296, 502 

D'Antona, F. & Mazzitelli, I. 1994, ApJS, 90, 467 

Desch, S. J. 2007, ApJ, 671, 878 



Dodson-Robinson, S. E., Veras, D., Ford, E. B. & Beich- 
man, C. A. 2009, ApJ, 707, 79 

Dullemond, C. P., Hollenbach, D., Kamp, I. & D'Alessio, 
P. 2007, Protostars and Planets V, ed. Reipurth, B., 
Jewitt, D. & Keil K. (Univ. of Arizona Press, Tucson), 
555 

Dvorak, R., Froeschle, C. & Froeschle, Ch. 1989 A&A, 
226, 335 

Eggleton, P. P. 2006, CUP, Evolutionary Processes in Bi- 
nary and Multiple Stars 

Figer, D. F., Najarro, F. & Kudritzki, R. P. 2004, ApJ, 
610, 109 

Figer, D. F. 2005, Natur., 434, 192 

Font A. S., McCarthy I. G., Johnstone D., Ballantyne D. 
R., 2004, ApJ, 607, 890 

Frank J., King, A. R. & Raine D. J. 2002, CUP, Accretion 
power in astrophysics 

Gammie, C. F. 2001, ApJ, 553, 174 

Geier, S., Edelmann, H., Heber, U., & Morales-Rueda, L. 

2009, ApJ, 702, L96 

Goodman, J. & Tan, J. C. 2004, ApJ, 608, 108 

Hayashi, C. 1981, Prog. Theor. Phys. Suppl., 70, 35 

Hosokawa, T. & Kazuyuki O. 2009, 691, 823 

Hollenbach, D., Johnstone, D., Lizano, S. & Shu, F. 1994, 
ApJ, 428, 654 

Isella, A., Carpenter, J. M. & Sargent, A. I. 2009, ApJ, 
701, 260 

Kashi A. & Soker, N. 2010, preprint, (|arXiv:0912.1439fr 2) 

Kalas, P. G., James R. C, Eugene, F., Michael P., 
Clampin, M., Kite, E. S., Stapelfeldt, K., Marois, C. 
& Krist, J. 2008, Sci, 322, 1345 

Kenyon, S. J., Bromley, B. C. 2009, ApJL, 690, 140 

Kratter, K. M. & Matzner C. D. 2006, MNRAS, 373, 1563 

Kratter, K. M., Murray-Clay, R. A. & Youdin, A. N. 

2010, ApJ, 710, 1375 

Krumholz, M. R., Klein, R. I., McKee, C. F., Offner, S. 
S. R. & Cunningham, A. J. 2009 Science, 323, 754 

Kuchner, M. J. 2004, ApJ, 612, 1147 

Kuiper, R., Klahr, H., Dullemond, C, Kley, W. & Hen- 
ning, Thomas 2010, A&A, 511A, 81 

Maeder, A., & Meynet, G. 1994, A&A, 287, 803 

Mayer, L., Quinn, T., Wadsley, J., & Stadel, J. 2002 Sci., 
298, 1756 

Mordasini, C, Alibert, Y., Benz, W. & Naef, D. 2009, 
A&A, 501, 1161 

Mullally, F., Winget, D. E., De Gennaro, S., Jeffery, E., 
Thompson, S. E., Chandler, D., Kepler, S. O. 2008, 
ASPC, 398, 163 

Musielakl, Z. E., Cuntz, M., Marshall, E. A., & T. Stuit, 
D. 2005 A&A 434, 355 

Najarro, F., 2005, ASPC, 332 58 



Nero, D. & Bjorkman, J. E. 2009, ApJL, 702, 163 

Perets, H. B. 2010 l|arXiv:1001.0581|) 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, 
J. J., Podolak, M., & Greenzweig, Y. 1996, Icar., 124, 
63 

Qian, S.-B., Liao, W.-P, Zhu, L.-Y., Dai, Z.-B., Liu, L., 
He, J.-J., Zhao, E.-G. & Li, L.-J. 2010, MNRAS, 401L, 
34 

Rice, E. L., Barman, T., Mclean, I. S.; Prato, L. & Kirk- 
patrick, J. D. 2010, ApJS, 186, 63 

Rice, W. K. M., Armitage, P. J., Bonnell, I. A., Bate, M. 
R., Jeffers, S. V. & Vine, S. G. 2003, MNRAS, 346L, 
36 

Schaerer, D. & de Koter, A. 1997, A&A 322, 598 

Schneider, J. 2010, The Extrasolar Planets Encyclopaedia 
( http: / / exoplanet.eu ) . 

Silvotti, R. et al. 2007, Nature, 449, 189 

Smith, N. & Owocki, S. 2006, ApJL, 645, 45 

Smith, N., Whitney, B. A., Conti, P. S., de Pree, C. G. 
& Jackson, J. M. 2009, MNRAS, 399, 952 

Shakura N. I. & Sunyaev R. A. 1973, A&A, 24, 337 

Stamatellos D. & Whitworth A. P. 2008, A&A, 480, 879 

Szebehely, V. 1980, CeMec, 22, 7 

Toomre A. 1964, ApJ, 139, 1217 

Vaidya, B., Fendt, C. & Beuther, H. 2009, ApJ, 702, 567 

Vlemmings, W. H. T., Surcis, G., Torstensson, K. J. E. 
& van Langevelde, H. J. 2010, MNRAS, 404, 134 

Weidenschilling, S. J. 1977, Ap&SS, 51, 153 

Wolszczan, A. & Frail, D. 1992, Nature 355, 145 



This 2-column preprint was prepared with the AAS IAT^X 
macros v5.2. 



